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Synergistic integration of electromagnetic (EM) and mechanical properties of metamaterials, a concept 
known as smart metamaterials, promises new applications across the spectrum, from flexible waveguides to 
shape-conforming cloaks. These applications became possible thanks to smart transformation optics (STO), 
a design methodology that utilizes coordinate transformations to control both EM wave propagation and 
mechanical deformation of the device. Here, we demonstrate several STO devices based on extremely auxetic 
(Poisson ratio — 1) elasto-electromagnetic metamaterials, both of which exhibit enormous flexibility and 
sustain efficient operation upon a wide range of deformations. Spatial maps of microwave electric fields 
across these devices confirm our ability to deform carpet cloaks, bent waveguides, and potentially other 
quasi-conformal TO-based devices operating at 7 ~ 8 GHz. These devices are each fabricated from a single 
sheet of initially uniform (double-periodic) square-lattice metamaterial, which acquires the necessary 
distribution of effective permittivity entirely from the mechanical deformation of its boundary. By 
integrating transformation optics and continuum mechanics theory, we provide analytical derivations for 
the design of STO devices. Additionally, we clarify an important point relating to two-dimensional STO 
devices: the difference between plane stress and plane strain assumptions, which lead to elastic 
metamaterials with Poisson ratio —1 and — °°, respectively. 



Transformation optics 1 (TO) provided a new way to design sophisticated electromagnetic devices using the 
form invariance of Maxwell's equations under coordinate transformations 1,2 . To implement the intricate 
gradient permittivity and/or permeability distributions designed by the TO theory, many researchers 
employed the novel concept of artificially structured metamaterials, which consist of deeply subwavelength unit 
cells coupling with electromagnetic waves collectively as an effective medium. From the first introduction of TO 
theory 1 , various TO applications have been suggested using metamaterials, such as cloaking 1317 , arbitrary light 
guiding 18,19 , extreme imaging lenses 20 23 , and other interesting approaches to manipulating light 24,25 . Amongst 
those applications, invisibility cloaks 1,2 contributed most significantly to the continued public curiosity and led to 
the expansion of related research fields 26,27 . The omnidirectional TO cloak 1 effectively compresses an object to 
appear as a singular point, making it invisible from all directions; however it requires a wide range of refractive 
index between zero and unity, which is very difficult to achieve in realistic media without significant loss or 
dispersion 3 5 . This grand challenge has been tackled by only a handful of experimental efforts in the microwave 6,7 
and optical 8 frequency domains. On the other hand, the carpet cloak device compresses an object to a flat sheet, 
which requires only modest ranges of material properties 9 . Moreover, for the TE polarization (out of plane E- 
field), a carpet cloak transformation can be implemented entirely without magnetic properties; consequently, 
they have been designed with large operational bandwidth and negligible attenuation 1017 . The absence of mag- 
netic response in carpet cloak is best understood using the conformal 2 and quasi-conformal 9 TO theory; it turns 
out to be closely related to the local isotropy of the underlying coordinate transformation 28 . Since exact conformal 
maps are very restrictive, the majority of devices based on locally isotropic metamaterials resort to quasi- 
conformal maps (QCM) 9 . QCM-based TO has already enabled numerous electromagnetic devices made of 
dielectric materials only, such as carpet cloaks 1017 , arbitrary-shape waveguides 18,19 , or lenses 20 , all of which can 
operate with low loss and broad bandwidth. 

If a quasi-conformal transformation optics (QCTO)-based device has a deviation of its boundary shape, for 
example, resulting from an elastic deformation, the device requires a complete redesign in order to preserve its 
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functionality. To prevent degraded performance and maintain original 
function under an external load or other elastic deformation, self- 
adjustable metamaterials are necessary. Recentiy, a concept of smart 
metamaterials was introduced, which was conceived to enable shape 
reconfigurability of QCTO-based devices 28 . It was suggested that, in 
two dimensions - that is, for devices having long extruded shapes and 
limited to in-plane wave propagation - so-called hyperauxetic (having 
Poisson ratio less than — 1) materials would be necessary for an exact 
implementation of this idea. Here, we present a rigorous derivation of 
Smart Transformation Optics (STO) in two dimensions, and clarify 
the effect of choosing a plane stress vs plane strain configuration on 
the mechanical properties of smart metamaterials. 

In a smart TO device, the boundary mechanical load should deform 
each unit cell so that the new distribution of electromagnetic prop- 
erties would implement the desired TO map; in our case, TO maps 
must be limited to conformal or quasi-conformal transformations. In 
general, it is very difficult to integrate TO and solid mechanics, 
because the stress and strain distributions of a loaded structure in 
mechanical deformation are solutions of complicated equations - 
the generalized Hooke's law, which is the governing equation of elas- 
todynamics for linear elastic materials 29 . By noting the fact that the 
sum of angle changes (a + p, in Fig. la) between unit cell's four sides 
is represented as shear strain (Y xy = a + fi) in solid-mechanical 
deformation, the quasi-conformality (a + P = 0) was numerically 
tested through multiphysics simulation (COMSOL) for 2D materials 
(in xy-plane) without longitudinal (z- direction) deformation 28 . Thus, 
the auxetic material with large absolute value of negative Poisson's 
ratio (PR, v), which is extremely difficult to realize, is the necessary 
building block of STO devices based on (quasi-)conformal maps: it 
provides nearly zero shear strain, causing the unit cells to compress or 
expand virtually isotropically regardless of the deformation. 



For the derivation of smart transformation optics, it is convenient 
to use explicit coordinates in the original flat space (X) and the 
mechanically deformed space (x = X+u), where the vector field u 
is the displacement field, as described in Fig. lb 29 . Displacement, or 
strain, enters the generalized Hooke's law, which determined the 
equilibrium configuration of the device given the boundary condi- 
tions and volumetric mechanical loads for the static deformations of 
linear elastic homogeneous isotropic materials, the complicated 
stress-strain relations reduce to the Cauchy-Navier equations of 
elasticity for displacements. The design of a transformation optical 
device is converted to the calculation of 2D Laplace's equations, 
V%x' = 0 in xy-plane, to evaluate the coordinate grids deformation 
during the transformation from a flat space (X ) to a distorted space 
(x 1 —X+ «) 3031 . This method suggests that the calculation of the 
material parameters for a TO device is equivalent to the computation 
of spatial deformation field governed by the 2D Laplace's equation 
with proper boundary conditions. Starting from the elastostatic gov- 
erning equation, after showing Laplace's equations of a deformed 
coordinate {x'), we analytically derive the general solutions of smart 
TO devices by integrating transformation optics and solid mechanics. 

The non-auxetic approximate smart cloak demonstrated in 
Ref.28, which was made of silicone rubber with a not-so-large dielec- 
tric constant (e^ = 2.88), required an additional structural part (a 
triangular patch) at the bottom of the carpet cloak to achieve the 
extremely small value of the transformation Jacobian, which was 
necessary to realize the desired effective permittivity range for cloak- 
ing 28 . The deformation range for smart cloaking, defined as the range 
of allowed elastic deformation while it can automatically implement 
a valid QCM, is also limited to a narrow range. However, to achieve a 
versatile - and exact - smart TO device, such as, for example, a fully 







Figure 1 | Schematic and material property change induced by elastic deformation, (a) Schematic of the elastic deformation of a unit cell, (b) The 
coordinate transformation between the coordinates before (X) and after (x! = X + 2) the elastic deformation, (c-d) The effective permittivity ratio 



generated by an elastic deformation versus the deformed area change (J = det(F) ), in comparison with the quasi-conformal-mapped permittivity 



ratio (l/J). (c) When the dielectric constant changes as£<j = 3, 10, 50 for a fixed volume fraction of fa = 0.5. (d) When volume fraction changes as =0.1, 
0.3, 0.5 for a fixed dielectric constant of £<j = 10. 
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flexible waveguide, it is not possible to use the patch method because 
deformations may appear at various, unpredictable different loca- 
tions in the waveguide. In addition, we would have to increase the 
deformable range of smart TO at the same time. Although it is known 
that a highly negative PR structure (v ~ —10) with extremely high 
permittivity is a good candidate for smart metamaterials, it is tech- 
nically very difficult to realize. Because high permittivity material is 
not easily available in nature, it is important to design an efficient 
auxetic structure that allows the highest effective permittivity out of a 
given dielectric material. 

Here, we demonstrate a versatile smart TO device by making any 
range of elastic deformation automatically implement QCM without 
any additional structural part. We fabricated an auxetic two-dimen- 
sional (2D) plate of v ~ — 1 which enables us to achieve the max- 
imized possible effective permittivity (e 1 ^ = j out of an elasto- 

electromagnetic metamaterial made of a specific material e d = 10.2 
(PTFE ceramic, ROGERS Corp.). The extreme auxetic property (v ~ 
— 1) and the compressible structure without empty space allow us to 
create an exact, wide-deformation-range QCTO device from a single 
sheet of an elasto-electromagnetic metamaterial, without resorting to 
the patching method used in Ref. 28. Our measured 2D E-field map- 
ping data at 7 ~ 8 GHz show that our auxetic structure has smart TO 
electromagnetic properties for elastic deformation, such as smart 
cloaking and smart arbitrary waveguiding. Using the elastostatic 
governing (Cauchy-Navier) equations in solid mechanics, we show 
that the deformed coordinates (x' =X+u) of a 2D plate, made by 
arbitrary elastic deformation, satisfy 2D Laplace's equations, 

V^x' = 0, for specific materials of v 1 or v 03 in plane stress 

or plane strain conditions, respectively. By integrating transforma- 
tion optics and solid mechanics, we analytically derive the general 
solutions of smart TO devices. 



8 2 u 2 



8 2 u 2 
sx t 2 + 8X 2 2 



fi- 



ll- 



In the limiting case of i 



-1, 



a 2 Ml 

dXi8X 2 
8 2 ui 



u 2 



(3) 



8 2 u 2 8 2 u 2 
8X 2 + 8X 2 2 



0, 



8 2 u x 

8X 2 ~*~ dX, 

i.e., the elastic deformation always satisfies 2D Laplace's equations in 
xy- plane, so that V^w = S/ 2 X (X + u) = V^x' = 0. Because we use sliding 
boundary condition (roller boundary condition in solid mechanics) 
in this case, the Laplace's equation along with this proper boundary 
condition enables us to consider arbitrary elastic deformations equi- 
valent to the quasi-conformal transformations 30,31 . By integrating TO 
and solid mechanics, this directly certifies that the elastic deforma- 
tion (inside xy-plane) of a plate with finite thickness automatically 
implements quasi-conformal transformation if the material is auxe- 
tic as v — > —1 when the longitudinal deformation (z-direction) is 
allowed. 

If we consider another 2D plate in xy-plane constrained to have no 
longitudinal (z-direction) deformation, it is represented by the plane 



strain condition in solid mechanics 
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It typically represents a cross section that cuts an infinite or very long 
depth such that we can ignore any end effects, in other words, a unit- 
depth model 32 . From Eq. (1), we get following relations as derived in 
Methods, 
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Results 

General solutions found so that elastic deformations automa- 
tically implement quasi-conformal transformations. The govern- 
ing equations of elastodynamic problems of linear elastic materials 
are the generalized Hooke's law, a it = C^fyj k, I = 1, 2, 3), that 
relates stress (a,j) and strain (e^) with the elastic modulus tensor 
(Cyj-;)- To represent the coordinate transformation made by an 
elastic deformation from an original flat space X ; to the trans- 
formed space x' j (=X i + Uj), we have to simplify the generalized 
Hooke's law by solving directly for the displacements (u t ), instead 
of stress and strain. For the elastostatic deformations of a linear 
elastic homogenous isotropic solid with Poisson's ratio (v) at 
uniform temperature (T = const), the governing equations reduce 
to the Cauchy-Navier equations of elasticity as following 29 , 



8 2 u k 



8 2 U; 
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where C^a is the elastic modulus tensor, u ; displacement vector, and 
p 0 mass density. 

If we consider a two-dimensional (2D in xy- plane) plate with finite 
thickness, the longitudinal (z-direction) deformation is allowed and 
it is represented as the well-known plane stress condition in solid 



The mechanical loads 



mechanics ( 033 = <t 2 i = ff i3 = 0, -z — = 0 

are assumed to be constant throughout the thickness and the thick- 
ness of the plate can vary 32 . From Eq. (1), we obtain following two 
equations as explained in detail in Methods, 
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In the limiting case of v — » — °°, the elastic deformation satisfies 2D 
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Laplace's equations in xy-plane, so that V x « = V X (X- 



= 0. This shows that the elastic deformation (inside xy-plane) of a 
plate automatically implements quasi-conformal transformation if 
the material is auxetic as v — * —=0 when the longitudinal deformation 
(z-direction) is not allowed. 

General solutions consistently produce zero shear strains. If elastic 
deformation satisfy quasi-conformal transformation, shear strain 
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is supposed to be 0 28 , which makes shear modulus 

°yJ)J 

(G = shear stress/shear strain) become very large compared to bulk 
modulus (K). In this limiting case of G2>K for linear homogeneous 
isotropic materials, Young's modulus (E) and Poisson's Ratio (v) 

1. Hence, the QCM 



result in E = 
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condition leads to the requirement of auxetic materials with v 1, 

if we allow the longitudinal deformation, i.e., plane stress condition. 
In our previous work, on the other hand, we considered a 2D system 
with longitudinal independence, constrained to the plane strain 



condition. For a material of given PR (v), the effective PR 
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under 2D plane strain constraint is derived to be v| 

explained in Methods. This constraint to a 2D system without 
longitudinal deformation requires highly negative v such as — 10 to 

, D v 
set v. a = close to —1. 

5 *ff 1 - v 



Design of optimized electromagnetic properties for auxetic TO 
devices. An auxetic metamaterial can be designed by employing an 
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elastic and electromagnetic crystal, made of an incompressible 
dielectric material (e^) with volume fraction of and air or free 
space (e a ~ 1). The elastic deformation can be defined by the 
deformation gradient tensor, Fjj = 8x! i /8Xj, and the effective 

permittivity ratio of this crystal after (s'fffj and before (?, re f) the defor- 
mation leads to e' eff / e ref = {e a + (e d -e a )/ d /det(F)}/{e a 

e a )/ d } 28 . When we change the dielectric constant as 6 d 
for a fixed volume fraction of/^ 

plotted the effective permittivity ratio 



3, 10, 50 
0.5, as presented in Fig. lc, we 



'■rcf 



versus the deformed 



area change (J = det(E)), in comparison with the quasi-conformal- 
mapped permittivity ratio (1/7). This implies that the higher dielectric 
constant matches more precisely to the TO rule. On the other hand, 
while we change the volume fraction as = 0.1,0.3,0.5 for a fixed 
dielectric constant of e d = 10, as given in Fig. Id, we also compared the 



effective permittivity ratio 



'■rcf 



with the quasi-conformal-mapped 



permittivity ratio (1//). This indicates that the higher volume fraction 
matches more precisely to the TO rule. Finally the higher permittivity 
and volume fraction are more desirable to realize auxetic structures 
which have self-adjusting TO properties for elastic deformations. 

Fabrication of an auxetic structure (v 1) allowing the highest 

permittivity for a given material. To fabricate a real auxetic smart 
metamaterial sample, we used a rigid dielectric material with high 
permittivity of = 10.2 (PTFE ceramic, ROGERS Corp.). We made 
an array structure of a unit cell with four vertical rigid rods, 
alternatively cross-linked at the top and bottom, as described in 
Fig. 2a. If we compress (stretch) one leg among the four legs along 
the plane, the top joint of four legs lifts up (lowers down). Thus, all 
the linked four legs are going to have the same tilting angles (0) at the 



side view (see Fig. 2b), resulting in a reduced (expanded) square at the 
top view in 2D xy-plane (see Fig. 2c). Because the lateral extension is 
the same as the axial extension, the cross-linked four legs can be 
considered as a unit cell of a 2D auxetic material of Poisson's ratio 
v ~ — 1 in xy-plane. To get enough range of effective permittivity 
variation for a TO device including a carpet cloak, we choose initial 
volume fraction asfj ~ 0.5 at the initial tilting angle of 0 o ~ 4° and 
the background permittivity of s re f ~ 5. If we elastically deform the 
sample, the tilting angle changes from 0 o to 0 and the effective 



permittivity of the deformed unit cell becomes 



1 



(e d -l) + l 



where Jacobian / = 



{h sin 9 + l 2 cos 0} 2 
{h sin Q 0 + l 2 cos 6 0 } 2 

li the vertical 



volume fraction 



length, and l 2 the 



{li sin Oo + h cos Oq} 2 
diagonal width of a rod. Figures 2d,e are the photographs of our 
200 mm X 100 mm sample with Zj = 10 mm, l 2 = 1.5 mm before 
and after an elastic deformation. We can clearly see the auxetic 
behavior of our sample, especially from the movie in 
Supplementary Information. Because we have a realistic bottleneck 
to find a feasible dielectric material of high permittivity, it is 
important that our auxetic structure can shrink into a dielectric 
structure without any empty space, which enables us to achieve the 

maximized possible effective permittivity ^s' e ^ = Saj made of a 

specific material with dielectric constant of e#. 

Experimental demonstration of a versatile smart TO device. To 

demonstrate a versatile TO device, our auxetic device is used in our 
experiments for both smart electromagnetic cloak and smart 
arbitrary waveguide at the same time. We place the sample in a 
microwave 2D E-field mapping apparatus and get the E-field 
mapping data of incident and scattered electromagnetic waves 





Z-l sinS + Z 2 cos 6 






Figure 2 | Schematic and photographs of an auxetic smart TO device, (a) Array schematic of a unit cell with four vertical rigid rods, alternatively 
cross-linked at the top and bottom. Right figure describes the changed shape after a shrinkage deformation induced by axial compression (green arrow in 
left figure), (b-c) (b) A side view and (c) a top view of a unit cell. I,, l 2 are the vertical length and the diagonal width of a rod, respectively, 
(d-e) Photographs of our 200 mm X 100 mm auxetic sample with Z t = 10 mm, l 2 = 1.5 mm (d) before and (e) after an elastic deformation. 
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through our sample. We deformed our auxetic sample by a metallic 
bump parameterized as y = h ■ cos 2 {xnl\QQ) in milimetres. 
Figures 3a-e are the electromagnetic cloaking results for 45° 
incident wave, while we change the bump height (h) as h = 0 mm, 
5 mm, 10 mm, 15 mm, and 20 mm, respectively. When we deform 
the auxetic structure by a bump, it redistributes the unit cells to 
satisfy quasi-conformality, consistently with our theoretical 
prediction (see the photographs of the first column). In the E-field 
mapping data at 7 GHz (see the second column) and 8 GHz (see the 
third column), the scattered waves simply reflect and do not split for 



a wide range of bump heights from h = 0 mm to 20 mm, in other 
words, the varying bumps are electromagnetically smart-cloaked. 
Without an auxetic cloak, as described in Fig. 3f, a perturbation of 
the metallic surface (h = 10 mm) causes a strongly scattered 
secondary beam with a power gap in a homogeneous medium of 
10 mm-thick silicone rubber plate (s d ~ 2.9) at 8 GHz. 

We also measured the electromagnetic wave propagation through 
the same auxetic sample, as presented in Figs. 4a-e, while we bend 
the sample at the angles of 0°, 30°, 45°, 60°, and 90° between the input 
and output facets (see the first column). The E-field mapping data 



7 GHz 8 GHz 




h - 10 mm 



Figure 3 | Experimentally measured E-field mapping of smart cloaking, (a-e) Experimentally measured electromagnetic cloaking results for 45° 
incident wave at 7 GHz and 8 GHz, while we change the bump height (h) as (a) h = 0 mm, (b) 5 mm, (c) 10 mm, (d) 15 mm, and (e) 20 mm, 
respectively, (f) Without an auxetic cloak, a beam incident on the bump of h = 10 mm in a homogenous medium of silicone rubber at 8 GHz. The 
photograph of our sample in each case is given in the first column. 



SCIENTIFIC REPORTS | 4 : 4084 | DOI: 1 0. 1 038/srep04084 



5 



7 GHz 8 GHz 




30° bending 




45° bending (8 GHz) 90° bending (8 GHz) 

Figure 4 | Experimentally measured E-field mapping of smart arbitrary waveguiding. (a-e) Experimentally measured electromagnetic wave 
propagation through the same auxetic sample at 7 GHz and 8 GHz, while we bend the sample at the angles of (a) 0", (b) 30°, (c) 45°, (d) 60°, and (e) 90 
between the input and output facets, (f-g) Without the auxetic sample, experimentally measured electromagnetic wave propagation through a 
homogenous waveguide made of silicone rubber at 8 GHz for (f) 45° and (g) 90° bending angles. The photograph in each case is given in the first 
column. 



are experimentally obtained at 7 GHz (see the second column) and 
8 GHz (see the third column). The phases of electromagnetic 
waves continuously maintain through the auxetic sample for any 
bending angles, i.e., "smart arbitrary waveguiding" behaviors. 
Figures 4f-g provide the electromagnetic wave propagation data 
through a homogenous waveguide made of silicone rubber medium 
at 8 GHz for 45° and 90 u bending angles. The internally reflected 
waves at the boundary interfaces of the waveguide interfere with the 
propagating waves, which produce phase fringe patterns of wave 
propagations inside the waveguide. These results experimentally 
show that, for an auxetic metamaterial made of a high permittivity 



dielectric material, the elastic deformation automatically implements 
quasi- conformality in the broadband microwave regime. 

Discussion 

We have made a 2D-auxetic metamaterial structure, whose unit cell 
consists of four vertical rigid rods, cross-linked at the top and bottom, 

and achieved Poisson's ratio v 1 for arbitrary deformations in the 

xy-plane. Since this structure exhibits smart elasto-electromagnetic 
behavior only when the dielectric inclusions have a very high dielec- 
tric constant, it is important to choose a geometry that minimizes 
empty space. Our choice of the high-index dielectric is the PTFE 
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ceramic with = 10.2 (Rogers Corp.). When we deform our auxetic 
carpet cloak sample with variable-height metallic bumps, our mea- 
sured E-field maps at 7 ~ 8 GHz indicate that the incident beam 
undergoes only a specular reflection without side lobes. This beha- 
vior is observed for a wide range of bump heights from h = 0 mm to 
20 mm. We also measured electromagnetic wave propagation 
through the smart waveguide sample, while we bend the sample at 
the angles of 0°, 30°, 45°, 60°, and 90° between the input and output 
facets. From the E-field mapping data at 7 ~ 8 GHz, the proper 
phases of electromagnetic waves are continuously maintained 
through the range of bending angles, i.e., the waveguide is entirely 
flexible. These experimental results show that, for an auxetic meta- 
material made of a high permittivity dielectric material, the elastic 
deformation automatically implements quasi-conformality in the 
broadband microwave regime. The higher permittivity and volume 
fraction are more desirable to realize auxetic structures which have 
self-adjusting TO properties for elastic deformations. 

In this work, using the elastostatic Cauchy-Navier equations 29 , we 
prove that the deformed coordinates of a 2D plate, made by arbitrary 
elastic deformations, satisfy Laplace's equations for specific materials 
of v = — 1 or v = — 00 in plane stress or plane strain conditions. The 
Laplace's equations along with sliding boundary condition allow us 
to prove analytically that the elastic deformation of a plate is auto- 
matically quasi-conformal. 

Methods 

Derivation of Laplace's equations from elastostatic governing equations. For plane 

d 

stress condition, (7 33 = c 31 = <r 23 — 0> — —0 an d stress-strain relation is 29 
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For plane strain condition, £ 33 
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Cauchy-Navier equations of motion, i.e., elastostatic equations become 
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If we insert all these coefficients into Cauchy-Navier equations, we get 
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In the limiting case of V — » — 1 , finally we obtain 2D Laplace's equation of u y as 
follows 
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From the stress-strain relations of plane stress configuration, we obtain each 
coefficient as followings, 
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E 



C25 — 0> Cl221 — 



2(l + v) 



^2221 — C1222 — C 2 e —C(a — 0, C2222 —£22 - 



E 

1^7 



02232 — ^24 — 0. 

After inserting all these coefficients into Cauchy-Navier equations, we get 

Ev 8?u\ E 8 2 u\ E 8 2 U2 E d 2 U2 



l-v 2 dx 1 8x 2 2(l+v)dxidx2 2(1 + v) dxf 1-v 2 8x 2 - 



= 0, 



d 2 u 2 8 2 u 2 
8xi 2 8x 2 2 



(1 + v) / d 2 u x 8 2 u 2 
(1 — v) \8x l 8x 2 8x 2 2 



(12) 



(8) 



d 2 Ui 8 2 ui d 2 Ui 
£-1111 — - +C2111 - — I-C1112 - — h 

OXi L CX\ OX2 OX\ cx 2 

d 2 Ul 8 2 u 2 d 2 u 2 

C-2U2 ^ — T +L-1121 « — T + W121 * « — + 

CX2 Z GJCi z CX1CX2 

5 2 U2 C 2 U2 5 2 U3 

W122 ^ 5 KC2122 7 + W131 - + 

t7X!C7X2 OX 2 ^ OXj 

5 2 U 3 C 2 W 3 S 2 «3 

W131 « a ^Cn 3 2^ i — +W1323 — T = 0. 

PX1DX2 OX 2 L 



For linear elastic materials, four rank (3X3X3X3) elastic modulus tensor 
C»kl is simplified as a matrix (6 X 6) C a ^ by a relationship of [if) — > a, (/c/) — > j9 
so that {(11),(22),(33),(23,32),(13,31),(12,21)}^(1,2,3,4,5,6), respect- 
ively 29 . Because the matrix C^p consists of components from the stress-strain 
relations of plane stress configuration in Eq. (6), we get each coefficient as 
followings, 



In the limiting case of v - 
follows, 



— 1, we also get 2D Laplace's equation of u 2 as 

d 2 u 2 d 2 u 2 



dx\- 



>0. 



(13) 



(2) Plane strain configuration: £ 3 3 = £ 3 i = £23 = 0, ^ — — 0. 

8x3 

For j — 1, we get coefficients of Cauchy-Navier equations, Cum 

Pq ~J^r = ^ an ^ er simplifying with - — — 0, as following, 

cx$ 



d 2 u k 
dxjdxj 



dt 2 

£1111 =£n — ■ 



/, . ^w, — ~^~\>^2U1 —Cmi — Cui2 —^61 —C\6 — 0>C]L3i — C]5 — 0, 

(l + v)(l-2v) ^ 

^2112 — £2121 —C^— , r, ^2131—^66—0, C\y 22 —C-]2 z 

2(1 + v) 

C 2 ]22 — Cq. —0, Cll32 — C\\ — 0, C2132 — C{A — 0. 



(l + v)(l-2v)' 
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After inserting above coefficients, the Cauchy-Navier equations finally become 
E(l-v) d 2 m E d 2 ui 



( 1 + v) ( 1 - 2v) ftt[ 2 2(1 + v) 8x 2 2 
E 8 2 u 2 Ev d 2 u 2 



2(1 + v) 8xi 8x 2 (1 + v) (1 -2v) Sxi 8x2 



= 0, 



5 2 wi d 2 m 
dx? dx} 



In the limiting case of v 
from Eq. (14). 



1 (d 2 U! d 2 u 2 
(l-2v) \dxj + dxidx 2 



(14) 



, D Ui 8 Ui 

c , we obtain Laplace s equation =- H T — >0 



dxjdxi 



For j — 2, we get coefficients of Cauchy-Navier equations, C j2 
p Q ^ — 0 after simplifying with ^— — 0, as following, 

£1211=^61=0, C2211 =£21 = 777 — \n ^TT' ^1212 =^1221 =^66 = T777 — \> 

(l + v)(l— 2v) 2 (!+ v ) 

^2212 = ^2221 — ^1222 — (?26 — C&2 — 0, ^1232 —^64 = 0, Cl231=C65=0, C2231 = 



2221 — ^1222 — 1-26 — (-62 



^26—^62 — 0, C2: 



-Cr. 



C 24 = 0. 

The Cauchy-Navier equations finally become 

Ev d 2 ui E 



(l + v)(l-2v)' 



d 2 ui 



(1 + v) (1 - 2v) 8x1 8x 2 2(1 + v) 8x1 dx 2 
E 8 2 u 2 E(l-v) d 2 u 2 



2(1 + v) dx}+ (l + v)(l-2v) dx 2 



= 0, 



8 u 2 d u 2 
~dx7 + l>xj 



C Ui 



8 2 u 2 

(1 — 2v) \dxidx 2 ' dx 2 



(15) 



, 8 U 2 d u 2 

, we obtain Laplace s equation _ + 2 — >0 

OX-i OXj 



In the limiting case of v - 
from Eq. (15). 

Derivation of effective 2D Poisson's ratio (v^fj in plane stress and plane strain 
configurations. Let's consider a planar material in xy-plane. To calculate 2D 
Poisson's ratio inside xy-plane, we compress the material in x-direction, i.e., g x (o u ) 
0, then calculate the planar Poisson's ratio defined as v^t = — ?, y j?, x . 



(i) 



Plane stress condition 

From the stress-strain relations, we get en 



£22 = 



(2) 



--Bx- E a*. - -, E 
Thus, the effective 2D Poisson's ratio becomes = — — — v. For the plane 
stress condition, the effective 2D PR is the same as the Poisson's ratio of the 
linear elastic homogeneous isotropic material. 
Plane strain condition 

From the stress-strain relations, we get En =£ x = — — — cr*, £22 = 
Ey= — V ^e"*" &x> Thus, the effective 2D Poisson's ratio becomes v 2 Jjj = 

By _ 



E x 



1-V 



For the plane strain condition, the effective 2D PR is different 

from the Poisson's ratio of the linear elastic homogeneous isotropic material. 
In this case, when v is highly negative number, such as — 10, is getting close 
to -1. 



Electric field measurements. Using a phase -sensing, near-field microwave scanning 
system within a planar waveguide, we map the TE field distribution of the scattering 
region, including the incident and scattered beams. The microwaves are launched into 
the planar waveguide from a standard X-band coax-to waveguide coupler. The sample 
was tested with 7-8 GHz waves, then capturing the TE field distribution in the xy-plane. 

Auxetic crystal fabrication. For a double- periodic and uniform auxetic crystal, we 
use PTFE ceramic (g^ — 10.2, ROGERS Corp.) with 1.2 mm side length square. The 
ceramic was cut into slices of 10 mm length, then alternatively bonded to each other 
at the top and bottom using scotch wood adhesive (AD6005, 3M Corp.). 
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